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SUMMARY 


In the process of analyzing the longitudinal frequency-response 
characteristics of aircraft, information on some of the methods of analy- 
sis has "been obtained by the Langley Laboratory of the National Advisory 
Committee for Aeronautics. In the investigation of these methods, the 
practical applications and limitations were stressed. 

In general, the methods considered may be classed as: (l) analysis 

of sinusoidal response, ( 2 ) analysis of transient response as to har- 
monic content through detennination of the Fourier integral by manual or 
machine methods, and ( 3 ) analysis of the transient through the xise of 
least-squares solutions of the coefficients of an assumed equation for 
either the transient time response or frequency response (sometimes 
referred to as cvirve-fitting methods) . 

The investigation has led to the following obser-vations: The cvirve- 

fitting methods (Donegan-Pearson and exponential-approximation methods) 
appear to be less critical to inputs having regions of low harmonic con- 
tent than Fourier methods and present the frequency response as analyti- 
cal expressions (transfer functions) . Fcrurier methods indicate charac- 
teristics of frequency response that may be missed in curve-fitting 
methods because of the limitations on the assumed form of the equations 
used in the curve-fitting methods. For manual calculations, the Ddnegan- 
Pearson method appears to be best suited for highly damped systems !in 
response to arbitrary control inputs, the exponential-approximation 
method appears "bo be best suited for lightly damped systems in response 
to step or short-pulse control inputs, and the Fourier method off era 
comparable results but requires lengthy calculations . Special machines 
for performing the Fourier analysis, such as the Coradl harmonic analyzer 
and the Foiirier sytheslzer, reduce the time required for the solution 
but do not offer particular improvement in accuracy over the us\ial 'man'ual 
methods . The use of punch-card calculating machines for the e-valudtion 
of the Fourier integrals appears to offer possibilities of more accurate 
resialts with a large reduction in time over the usual manual methods . 
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INTRODUCTION 


In recent years, a large number of methods have "been advanced for 
the purpose of obtaining frequency-response data, transfer functions, 
and stability coefficients from flight tests by using control inputs of 
arbitrary shape. As pointed out in reference 1, the data obtained from 
the application of these methods axe of great value to aircraft and auto 
pilot designers as well as to designers of other electronic airborne 
equipment for the combination of their individual products into a stable 
working vinit . 

A resum^ of methods and progress to date wit^ reference to dynamic 
flight testing is presented by Milliken in reference 2. Although all 
these methods appear to have advantages and limitations, some methods 
have gained popularity with various groups whereas other methods remain 
comparatively tmused. A number of methods have been examined and used 
by the Flight Research Division of the Langley Aeronautical Laboratory 
in an effort to determine which methods to adopt~^n establishing the 
transfer functions of the various aircraftncmdergoing dynamic flight 
tests . This program has offered a certain amount of practical experi- 
ence in the use and limitations of the methods, and it is believed that 
this experience may be of value to others engaged in obtaining the fre- 
quency responses of .aircraft. No attempt has been made, to examine all 
the known methods of analyzing dynamic responses and omission of any 
method is not intended to Imply lack of merit. 

A brief review of the methods examined is offered in this paper, 
together with references to their derivations and examples of their 
application. Three .types of aircraft, a fighter, a transport, and a 
free-fall model, were used for these examples. The examples are con- 
cerned with the short-period longitudinal mode of the airplane, which 
is usually a well-damped mode defined by a fairly simple transfer func- 
tion. Thus, the comparisons of methods presented herein are made solely 
on the basis of resTilts obtained from analysis of this longitudinal mode 
It is recognized that complicated, oscillating systems may be analyzed 
and greater accuracies obtained by all the methods reported at the 
expense of a more extensive, analysis . The comparison of methods pre- 
sented herein, however, may be altered when applied to more complicated 
systems . 

The methods are discussed with regard to the time required, the 
means for facilitating their use, and the limitations on their applica- 
tion. Some opinions presented are not directly substantiated by quanti- 
tative results but are based on experience in the use of-the methods. 

The results obtained are compared to give some indication of the rela- 
tive accuracy of the methods, exclusive of any inaccuracy in the 
measurements . 
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SYMBOLS 



3 = F 

Ki,Ka 

k 

M 

m 

n 

P 

P 

Q 

<1 

T 


coefficients of longitudinal transfer functions of an 
airplane (composed of stability parameters) 

amplitude 

coefficients of cuLic approximation to transient 

damping e^qjonent in e®’'^ 

dandling coefficient, -2a 

arbitrary coefficient 

differential operator, d/dt 

natural logarithmic base, 2 . 7 I 828 

pressure altitude, ft 

coefficients of eq-uation relating airplane longitudinal 
response to step-control motion 

real and imaginary terms of Fourier integral, respectively 

stiffness coefficient, a^ + 

Mach number 
numerical Integer 

numerical integer or normal acceleration, g units 
period, sec 

Laplace transform variable 
quantity as a function of frequency 
quantity as a fiinction of time 
time variable of integration, sec 
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^ 1/2 

'^ 1/20 

t 

tp 

tr 

X 

y 

a 

7 

B 


e 

T 

0 

(U 

o)n 

O 

Subscripts ; 

e 

I 


time to dem^ to one-half amplitude, sec 
time to damp to one-twentieth an 5 )litude, sec 
time, sec 

time required for oscillation to reach initial peak 
from time of step-control input, sec 

time interval over which-ciarve analyzed 

abscissa of response curves 

ordinate of- response curves 

angle of attack, deg 

flight -path angle, deg 

control deflection, deg 

Danpixig coefficient 


damping ratio. 


Critical damping coefficient 


pitch angle, deg 
time lag, sec 
phase angle, deg 

displacement coefficients of Coradi harmonic analyzer 


frequency, radians /sec 

•undaniped natiiral frequency, radlans/sec 

forcing frequency, radians/sec 


elevator 

input 


0 


output 
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o 

magnitude 

at t = 0 

ss 

magnitude 

at steady state 

00 

magnitude 

at Infinity 

€ 

error 


0,1, 2, 5, . • .n 

numerical 

integers 


The absolute value of any term is denoted by | | . 

BASIC CONCEPTS 


It Is assumed herein that the reader Is familiar with the concepts 
and application of the Laplace transform to linear systems. A presen- 
tation of this method may be found in reference 3* 

The frequency response of a dynamic system defines its steady-state 
response under the influence of an input applied in the form of a sinu- 
soidal oscillation of constant amplitude and period. An analytical 
ejspression which defines the freq.uency response throu^out the frequency 
range is, when defined in terms of the Laplace transform variable p, 
the transfer function of the system. The transfer function not only 
expresses the freq.uency response but it may be said that, for linear 
systems, any arbitrary input fiinction operated on by the transfer func- 
tion determines the variation in the output function. Conversely, if 
the input and output are known, it should be possible to determine an 
analytical expression which relates the two, that is, the transfer f\mc- 
tion. The present paper is concerned with several methods of obtaining 
the transfer functions of aircraft from measiired Inputs and responses. 

The methods presented herein, in general, may be divided into two classes 
methods that first determine the frequency response of the system and 
methods that determine the transfer function without the determination 
of the frequency response. 

The MCA sign convention, as shown in figure 1, assumes elevator 
trailing edge down as positive. Therefore, a positive elevator deflec- 
tion will, in general, produce negative static responses. In order to 
conform with the usvial practice of plotting frequency-response data, 
phase angles have been shifted l80° (that is, zero phase angle at zero 
frequency) . 
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DESCRIPTION AND DISCUSSION OF METHODS 
Sinusoidal-Response Method 


Of the several possible ways to obtain the frequency response of 
a system, an obvious way is to oscillate sinusoidally a control surface 
at a constant airiplit\ade and frequency until a steady-state response of 
the aircraft has been obtained and measure the amplitude and phase rela- 
tionship between input and output' sine waves. The process may then be 
repeated throughou-fc-the frequency range of interest. The theoretical 
application of this method to the determination of the coefficients of 
the transfer function is given by Greenberg in reference k, and a graph- 
ical method of detenainlng transfer functions from frequency-response 
data is given in reference 5* 

The sinusoidal-response method requires the least computation time 
and the most flight time of the methods reported herein. In an effort 
to reduce the large amount of flight testing required, a number of 
simplifications from the usual technique have been attempted. One pro- 
cedure that was investigated involved obtaining sinusoidal-response 
data by contlnuoiisly recording the controlled -input , and the response of 
the airplane while slowly changing the frequency of the input to cover 
the range of- frequencies desired. Appendix A presents an estimation of 
the errors encountered at several values of- rate of change of forcing 
frequency when such a frequency-modulated input is applied to a dynamic 
system defined by a second-order lag. From this analysis and also from 
flight results, it appears that, for systems having near critical damping, 
satisfactory results may be obtained. In addition, it" appears that a 
human pilot may generate an adeqtiately near sinusoidal input without the 
aid of special equipment, particularly if he has a fairly precise indi- 
cation of the amplitude of his stick motion. A typical record obtained 
by using these techniques is presented in figure 2 . The deviation from 
a pure sinusoidal input is obvious although the filtering supplied by 
the airframe results in a nearly sinusoidal response. Jones and 
Sternfield in reference 6 outline a method for determining the ampli- 
tude of an equivalent sine wave when the actual periodic wave has an 
irregular form. In general, however, it has been possible to obtain 
results consistent with the accuracy of- the measurements by fairing the 
peaks of the oscillations in the input and output and obtaining the 
double amplitude of these quantities from the fairing by averaging over 
a nxmiber of successive half-cycles. The mean value about which the 
oscillations occur is established from the fairing of the peak ampli- 
tudes, and the time lag of the output behind the input is determined by 
averaging the lag read along this mean value over a number of successive 
half-cycles. The period of the oscillations is similarly obtained by 
averaging. The method of measurement of these quantities is Illustrated 
in fig\rre 2. 
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The in55ortaijce of averaging over successive half -cycles, vhen 
establishing the time lag. Is shown In figure 2 where, in some cases, 
the velocity of the input Is considerably different in one direction 
than in the other with the resiilt that the time lag read at one point 
will differ considerably from that read one-half period later. The 
average, however, appears, in most cases, to represent adequately the 
actual time lag. 

The frequency-response parameters may be determined from the 
averaged values through use of the following relations: 


Frequency = 


2it 

Period 


( 1 ) 


An^jlitude ratio = 


Amplitude of output 
Amplitude of input 


( 2 ) 


Phase angle = 


560 — 

Period 


(5) 


Flight records for a fighter airplane were obtained with the pilot 
manually applying an approximate sine wave of varying frequency to the 
elevator. A sample of the flight data is presented in figure 2 and 
some pertinent geometric characteristics of the airplane and the flight 
condition under which the data were obtained are listed in table I. 

Data points obtained by the foregoing analysis are presented in fig- 
Tore 3 together with a suggested fairing. The scatter is considered 
typical for this technique (pilot- induced input) . These data points 
represent portions of two flights of the fighter airplane and a recording 
time of about 25 O seconds. A typical example of the time required to 
reduce the flight data to a frequency-response curve by this method is 
shown in table II. typical times are also presented in the table for 
other methods to be discussed subsequently. 


Fourier Analysis of Transient Response 

Another well-known method of determining the frequency response is 
to dete rmin e the coefficients of the Fourier transform of the input and 
output functions over a frequency range by analyzing the response (as a 
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function of time) of the aircraft to an arbitrary input. The process 
is indicated hy the expression 


^0 

% 


(Joi) 



(4) 


which represents the ratio of the Fourier integral of the output to the 
Fourier integral of the inpiib. The derivation and s"everal applications 
of this method are reported in references 4 and 7 . 

This method^ as well as the sinusoidal method previously described, 
gives data points in amplitude ratio and phase angle at a number of dis- 
crete frequencies through which a fairing can be made. Transfer func- 
tions may be obtained from the frequency-response curves by the same 
methods as were mentioned for the sinusoidal response. 

Integration of the Fourier integral offers a choice of methods 
which may be divided into two general categories: (_l) methods which 

divide the transient- into finite intervals, approxiiiiate the curves 
within each interval with an analytical expression, perform the indi- 
cated integration analytically, and sum the real and imaginary parts of 

these integrals; and ( 2 ) methods which express e'”'^^ in trigonometric 
form, multiply these sine and cosine functions by the value of q(t) 
at corresponding times, and integrate the product curves to determine 
the real and imaginary terms of the Fourier integral”. 

Solution of the Fourier integral by either approach involves the 
judicious choice of time Intervals. For the analytical representation 
method, a choice in the form of the analytical expression must also be 
made. As the chosen expression becomes more complex, the accinracy of 
the determination of the frequency response generally Increases but, 
as the work Involved likewise Increases, a compromise usually is 
necessary. 

The following are severa,l methods which have been studied and 
illustrate the various approaches to the Fourier tr^sfomiation. 

Manual Method: Analytical Integration within discrete Intervals o: 

cubic representations of a transient .- A method of representing a tran- 
sient for solution of the Fourier integral, as developed by Ordway B. 
Gates, Jr., of the Langley laboratory. Involves the division of the 
transient into discrete time Intervals chosen to facilitate acc\irate 
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approximation of each portion of the transient hy cubic (or lower-order) 
polynomials. The Fourier integral will then be 


Q(oi) 



, . -Jcut 
q.(t)e dt 


- 1 ’o'*' 





q.^(t)e“*^^dt + . 





q.j^(t)e "^^t 





q^g(t)e *^^dt 


(5) 


where 


(ln(t) = a^t^ + b^t^ + c^t + dn 


The values of the coefficients a^, b^, Cj^, and d^^ for any 

given n may be determined from the characteristics of the transient 
within the interval t^ to For the general case, the interval 

is subdivided into thirds and values of the transient q.(t) at these 
dividing points afford fo\jr cubic equations having four unknowns 
(an^ bn, and d^) . The advantage of using equal divisions within 


the interval tjj to is the ease of the solution of the four 


equations by means of "successive subtraction." (See illustrated 
example in appendix B.) If, however, the slope of the transient is 


zero 



within the inteirval. 


this condition should be used as 


well as the value of the transient q(t) at this point. The coeffi- 
cients thus determined give an equation that may be used to check the 
fit of the transient by the expression before further work is Initiated. 


This approach to the evaluation of the Fourier integral may be 
expressed analytically as follows: 


n=3S 

Q(oj) = ^ Qn(oj) = K^(co) - JK2(cs) 
n=0 


(6) 
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where 


^rM-1 




t‘^ + c_t + d„)e ‘^^dt 


“n 


and 


"^^n+l 


*n 


t^e'^'^^dt = - + 6jtoj - - jt^aj^)e"'^^ 


‘'n 


OJ 


Wl 


^n f t^e"‘^‘*^dt = ^(2 + 2jtcD - t^oj^)e“‘5'‘^ 

>5 ^ 


^n+1 




r^wl -Jolt Cn , . -Jta> 

Cn / te dt = — (1 + Jto))e 

'^tn <x>^ 


^n+1 


“n 


r*»i . is. e-l’” 

< *“ 




“n 


The substitution of discrete values of frequency o) gives the real 
and imaginary terms of the Fourier transform of the time transient q(t)-, 
and the relationships of anplitude and phase are given as. 

A = 1 + Kg 


0 = tan-1 

Kl 

As an illustration of this method, a numerical example is presented 
in appendix B. A method of this type is not very adaptable to machine 
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methods because some discretion is required in the subdivision of the 
transient, and the time intervals are not necessarily equal. A small 
amount of trial and error may be required in the choice of these time 
intervals, particularly in the vicinity of points of Inflection. The 
cubic representation has the merit of being the lowest-order polynomial 
to contain a point of inflection. This approach to the solution of the 
Fourier integral has the advantage of providing an analytical represen- 
tation that may be directly compared with the transient and of providing 
an exact analytical integration. On the basis of comparable accuracy, 
this approach is in many Instances shorter than the classical numerical 
Integration method that follows . 

Manml Method; Numerical integration of product c\jrves of 
q(t)sin (iit and q(t)cos cut. - The us\aal manual method, which also is the 

basis of some of the machine methods, requires a large amount of graphical 
or nxmerical integration because no attempt is made to obtain a contin- 
uous analytic expression for the transient until it has reached a steady 
state. If the input is restricted to a simple analytical expression 
(for example, a step input will have a constant value from zero time to 
infinity) , the graphical or nimierical Integration of the input is elimi- 
nated and the time required for the solution of equation (ii-) is roughly 
reduced by one-half. For the purpose of graphical or numerical integra- 
tion, equation (it) may be reduced to an expression involving real inte- 
grals of the form 




q(t)cos cut dt - 


j 



q(t)sln cot dt 


= K^(co) - JK2 (co) 


(8) 


where the frequency-response relationships are given by equation (7) . 

The numerical and graphical method for the solution of the Fo\jrier 
integral and an example of its use is given in reference it by Greenberg 
and a more cxamplete discussion is given by Schetzer in reference 8. 

A rule of thimib for choosing the proper time Interval in analyzing 
flight data has been suggested by experience gained in the use of this 
technique. The rule is restricted to the methods of integration adapt- 
able to the manual methods, for example, Simpson's three-point rule. 

In general, a chosen time interval At will give reasonable results 

up to a frequency of so that a time interval of 0.10 second may 


5 At 
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be expected to produce good results up to a frequency of about 6 rad.ians 
per second. The proper choice of time intervals is obviously dependent 
on the character of the input and output; however^ the foregoing rule 
has been found ustially conservative exce^~ in cases of extremely erratic 
variations in the input and output. The highest frequency at which 
reasonable results might "be expected will also depend upon the choice 
of forms of integration, a superior integrating method affording lise of 
greater time intervals for conpaxable results. 

When Simpson's three-point rule for numerical integration or a 
planlmeter are used for this method, the time required may be estimated 
by another rule gained from experience. For a typical, case where the 
short-period longitudinal response to a step or pulse input is analyzed, 
the time required to obtain the amplitude and phase angle of the output 
(one) function at 8 frequencies by ^lsing 24 data po^ts has been found 
to be about 6 man-hours for an experienced user. Tabulated values of 
sin cut and cos cut were used and the time required to make these tabu- 
lations was not included in the estimate. 

Punch- card method (IBM) .- Certain International Business Machines 
(described subsequently and referred to as IBM machines) offer a time- 
saving solution to the process outlined in the previous section with 
usutally more accurate resiilts over a greater frequency range since use 
of more complicated and accurate methods of mmierlcal integration are 
feasible. Weddle's seven -point rule (ref. 9 ^ P* 12 $) as derived from 
the Newton-Cotes quadrat-ure formula has been employed and is an example 
of such a method. 

By using essentially the same procedure as the msuiual method which 
integrates the product crurves q(t)sin cut and q(t)cos cut, a set of 
"master" cards are punched which define the values of the cosine and 
sine functions for the values of cut— selected and also define the 
numerical integration process used. Since values of cut determine 
the values of the trigonometric f\mctlons piinched on the cards, the 
frequency range to be evaluated may be varied by changing the time 
interval- in inverse proportion. Cards are likewise punched for the 
time functions of input and output (and these mvist obviously be punched 
for each separate analyzed record) . 

The calculations involved in the Fourier analysis method as per- 
formed on the IBM machines that are available at the Langley Laboratory 
are as follows; — 

(1) Time response data are perforated onto IBM cards by using a 
card punch . 

(2) Correct transcription of data onto cards is checked by a 
varif ier . 
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(5) Original deck, of cards representing time-response data is repro- 
duced, one deck for each, frequency to he analyzed, hy using a reproducer. 

(i^-) Integrating factors and trigonometric functions are transferred 
from the "master" deck to each deck obtained from step (5) t>y using the 
reproducer . 

(5) Product functions (q(t)sin oit and q(t)cos cot) are obtained by 
an electronic calculator . 

(6) End corrections and integration corrections are applied by using 
a sorter and the electronic calculator. 

(7) Fourier summation of terms obtained for each frequency are made 
by using a tabulator (alphabetical accoimting machine) . 

(8) Summary cards of amplitude and phase relationships are obtained 
on the electronic calculator. 

(9) Final frequency data of input and output functions are typed out 
by the tabulator. 

Some aspects of this process as applied to lateral responses having steady- 
state oscillatory responses are described in more detail in reference 10. 

The time required for this IBM equipment to perform the operations 
Indicated has been fo\md to be 5 machine-hours for the determination of 
data at 12 frequencies for one f\mctlon by using one set of machines and 
214-1 data points ( 12-second records using 0.05-second intervals). This 
time was averaged over several performances and included all checks and 
correction of mistakes. 

Method using the electromechanical Fourier synthesizer .- The elec- 
tromechanical FoTurler synthesizer, originally built and used by the 
Massachusetts Institute of Technology to produce transient-response 
curves from frequency-response data, was designed to perform the fol- 
lowing operation (see ref, ll) : 

') = YL ^ cos(np +0n)- 

n=l,2,5,... n=l,2,5,... 


J 



sln^np + 


( 9 ) 


n=l,2,5,. 
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where 

aitqplitude of nth harmonic 
phase angle of nth haimonic 
P angular displacement of fundamental 

The application of the Fourier synthesizer to the evaluation of the 
Fourier integral may be seen from the following derivation. 

The general form of analysis assumes that any arbitrary curve of 
input or output may be represented by a series of step fxmctions with 
a constant finite time lag between the steps: 



to t^ t2 t^ ti^. t^ tg 


Time, t 


The step approximation is the same as that used in the analytical method 
of obtaining the frequency response from a time response to a step input 
as presented in reference 12 and also the same as the extension of this 
analytical method to an arbitrary input as presented in reference I 5 . 
Thus, the Fo\irler Integral of an arbitrary function in time q.(t) may 
be approximated in the form 
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Q(cu) = 


Jtn - 


2 


Ai2_(t)e *^^dt 


poo 


1 . At 

Ut2-^ 




e dt 


• + 


tn- 


At 


, . -Jojt 
Ai^(t)e dt 


( 10 ) 


Therefore , 



(11) 


In trigonometric form 


Q(ai) 


i 7 

03 Z_ 


A3,jj^(t)sin^na3 At - o) 



1 ^ Ai^(t)co8fm3 At - 03 

n=1^2^5^ • • • 


( 11 a) 


This relationship, as can he seen hy con:5)arison of eq-uation (lla) with 
equation (9), may he handled hy the Fourier synthesizer. 

The nijimher of points that may he used conveniently to represent the 
time-response curve is determined hy the number of resolvers available 
in the machine to simulate the convolution process. The machine inves- 
tigated, employing 2 h resolvers, required 14- to 8 hours to obtain the 
frequency response of a system from any arbitrary input and output that 
may he represented hy 2 h eq-ually spaced steps . The frequency data are 
presented hy the machine as cvirves of the Fourier transforms of the real 
p.nt^ imaginary coefficients. From these curves, values at any momher of 
frequencies may he chosen for the determination of phase angle and 
amplitude . 
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In view of the fact that the Fourier synthesizer utilizes 2h equally 
spaced steps, it is limited to transients that may be adeqiiately approxi- 
mated thereby. Although the frequency range plotted by the synthesizer 
is from 0 to n/At radians per second, the results do not appear to be 
accurate to any higher frequency range than is quoted for the mamaal 

method which uses Simpson's rule f(a = — — . 

\ 5 At/ 

Coradl harmonic analyzer .- The Coradl harmonic analyzer (ref. 14) 
is a semimanually operated tracing machine which by the use of several 
rolling spheres may be \ised to evaluate the Fourier integral of a func- 
tion. The model Investigated (Dent-Draper Model, Rolling Sphere Type, 
Mico Instrument” Company, Cambridge, Mass.) employed five spheres which, 
through use of various gears, may measure the harmonic content of a 
curve within a range of 1 to 50 harmonics. Details of the operation of 
an earlier model of the Coradl harmonic ana3^zer are presented in 
reference 15. 

Throiigh the use of the Coradl harmonic analyzer, the time transient 
is traced from the point of initiation (initial conditions zero) to the 
point of steady-state response and the following integrals are evaluated: 


'Ifl 



cos cut d 



(I2a) 


^2 = 



sin cut d 



(12b) 


These integrals are proportional to the real and imaginary part of the 
Fourier integral of the cvirve being analyzed 


Q(<u) = f q(t)e”'^^dt 

The proportionality may be seen by integrating the Fourier integral by 
parts to change variables so that 


Q(oj) 


4(t) -jo>t “ ^ 1 

e + — / 

Jtw 0 



( 13 ) 
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■where for all. practical purposes the first term is zero. In the use of 
the Coradi harmonic analyzer, the term q(t) is plotted along the y-axis 
therefore, equation (15) inay he written as follows: 


Q(q) = — r “ e”'^\[y(tT] 
J n 


Jq 


icaJo 


'^00 -jolt 
/ f 
^'^Ybs 


e -dy+i /■‘'“e-'“\[;(t3 


(li^) 


where the second tenn is zero since = Yoo" trigonometric, form, 

equation (li^) becomes 


1 r^BB 


/"‘y's g 

Q(oj) = - ^ J sin cut d[y(t^ ^ J (^ 5 ) 


so that substitution of eqmtions (l2) into eq-uation (15) gives 


-to 


’t'l 

cdC 


(16) 


where cjo = and t-r is defined as the time interval over which the 

tr 

curve was analyzed. The constant l/C is the scale factor between the 
function being analyzed and the displacement of the dials of the analyzer 
In general, because only the ratio of output to input is desired, the 
individual scale factors need not be con5)uted provided both quan'titles 
are plotted to the same scale. 

This analyzer appears to produce the Fourier coefficients within 
an accuracy dependent upon the kinematic accuracy of the machine (which 
is primarily affected by slipping of the rollers but also to some extent 
by wear) and the ability of the operator and machine to follow exactly 
the trace being analyzed. The operator is required to follow the curve 
in the direction of the ordinate q(t) while the machine, operated by 
a microswitch, automatically traverses along the abscissa t. Accxarate 
tracing becomes diffiaalt when steep slopes (large values of dq/dt) are 
experienced, and a certain amoimt of roller slipping and h-uman error 
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should be expected. This inaccuracy is alleviated somewhat hy averaging 
the values obtained from three or more repeated tracings . Experience 
with the machine has indicated that the accuracies obtained are about 
the same as those obtained by the manual methods . 

The time required to obtain the Fourier coefficients of one func- 
tion q(t) for 15 haraionlcs has been found to be about ij.. hours for an 
experienced operator. This estimate includes the time to aline correctly 
the axes of the curve with the machine, to connect the correct set of 
gears for each 5 harmonics, and to trace the curve three times for each 
set of 5 harmonics. 

A consideration in the use of this machine is its ability to produce 
the Fourier coefficients in a comparatively short time, particularly with 
erratic functions that" would require very smalL time intervals to repre- 
sent accurately the function for use by other methods. A point worth 
noting is that, for erratic functions, the average of several tracings 
should produce a more reliable resialt. In the use of the Coradi harmonic 
analyzer, the limitation that the function must reach steady state still 
applies . 


Curve-Fitting ffethods 

In the methods herein called curve-fitting methods, the form of the 
treinsfer f\mction is directly or Indirectly assumed and the coefficients 
of the transfer function are determined by least-squares methods or a 
combination of least-squares and direct-conputabion methods . With a 
number of these methods, the analytical expression called the transfer 
function is obtained without first obtaining frequency-response data. 

The Donegan -Pearson method requires a direct-assunptlon as to the 
form of the transfer function and solves for the coefficients by sub- 
stituting into the transfer function the input and the output time 
functions and their integrals. On the other hand, the exponential- 
approximation method solves for the coefficients of analytical expres- 
sions which approximate the time histories of the input and output 
functions. The transfer function is then established by taking the 
Laplace transfonn of these analytical functions. 

In either the Donegan-Pearson or the exponential-approximation 
methods, the order of the expressions used to approximate""either the- 
transfer function or the input-output time histories is imlimited. 
Therefore, the requirement that the form of the transfer function be 
assumed would appear not to be particularly restrictive other than that 
the system be linear. In practice, however, the coirputation Involved 
in the least-squares procedure increases rapidly with increase in the 
order of the equations and the eqmtions tend to become progressively 
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more ill-conditioned. The general practice therefore has heen to assume 
a form for a given transfer function, as would he predicted from the 
stability theory, and this practice, in general, neglects low-freq.uency 
(phugoid) modes and possible high-frequency modes due to structtiral 
elasticity in order to hold the order of the equations to a minimum. 

Such proced-ures do not afford detection of these modes from flight-test 
data in cases where these modes are Important unless a form of the trans- 
fer function is assumed in advance which include such modes. In contrast 
Fourier analysis will detect all details of the frequency response which 
are within the accuracy of the measxirements and the calculation procedure 
The forms of the longitudinal transfer functions usually assumed in con- 
junction with the curve-fitting methods are: 

S Ap2 + Bp + C 


D0 Gp + H 

® Ap^ + Bp + C 

n _ Lp^ + + N 

5 Ap2 + Bp + C 

where the substitution of Jo) for the Laplace transform operator p 
gives the frequency response of the system. 

Exponential-approximation methods .- Since the response of a linear 
system to a step or impulse is a sum of exponentials, an obvious method 
for fitting airplane time responses is the choice of exponential terms. 
The nimiber of exponentials is selected so that the Laplace transforma- 
tion will give the same polynomial e3q>ressions as obtained from stability 
theory. Although this method can be applied to any input that has a 
Laplace transform, it is most suited for application to responses to an 
approximate Impulse, a step, and an approximate step input. The practi- 
cal "approximate step," as compared to the theoretical perfect step, may 
have a small but finite lag in reaching steady state and may have a small 
undershoot or overshoot. (See ref. l6.) The response eqioation to a step 
may be represented for the case of the short-period longitudinal mode of 
motion of an airplane by the form 


(it) 


q(t) = q^ + sin cot + Jq cos cut) 


(18) 
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where a is the damping exponential^ as is the freq.uency^ and J2 

are the coefficients of the in-phase and out-of-phase conponents of the 
response, and q.(t) is a function of time that expresses the response 
of the system. 

In the Laplace, transformation of the general form given hy equa- 
tion (18), the denominator of the transfer function would he given hy 


jp - (a +-Ja>^[p - (a - jos^ = p^ + hp + k 


(19) 


For lightly damped systems where the period and time to danp to one-half 
amplitude may he read directly from the response records, it" has heen 
found that direct“"calculatlon affords an accurate and rapid means of 
obtaining h and k. In the use of this method it follows from equa- 
tions (18) and (19) that 


b = = _2a (20) 

T / 

^ 1/2 


O.J+8 59.48 

k a ■' + — 

where '^ 2./2 time required for the oscillation to damp to one- 

half amplitude and P is the measured period of^the oscillation. 

This process, where the coefficients of the transfer function are 
computed from direct measurements of the flight records, has heen used 
extensively for the case of rocket and free-fall test models since these 
test models, in general, exhibit the low damping and high natixral fre- 
quency which enable this approach. 

Once a and (a are determined, values of J^_ and J2 appearing 

in eqviation (18) may he obtained from the time history. In Instances 
where the steady state is adeq\aately defined, direct conputatlon of J2 

is afforded. In the analysis of the response of an airplane in angle 
of attack and pitching velocity over short periods of time, the coeffi- 
cient J2 must he negatively eqtial to the value o^ the- response at 


+ (O 


( 21 ) 
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steady state. For the response in normal acceleration, however, a step 
input of the control surface causes an effective instantaneous change 
in load on the tail which, in turn, produces an instantaneous jun^ in 
the normal-acceleration response of the aircraft. An illustration of 
this type of response is shown in figure i4-. The relationship among 
the instantaneous change in acceleration at t = 0 , the steady-state 
acceleration, and J2 is shown hy eq.uation (18) for t = 0 so that 

a(o) = + Jg (22) 

With the coefficient J2 th\is established, a possibility for 
direct calc\iLation of may be indicated for the case of a step 

input and lightly damped systems by writing equation (18) in the form 


Q-(t) = 9 oo (25) 

and 

^ = tan“^ — = t (24) 

\2 P/ 

where P is the period of the osci 3 Ja,tion and tp is the time req.uired 

for the oscillation to reach a peak after initiation of the step input. 
The relationship among 0 , P, and tp may be visualized by reference 

to figure !<-. 

The method of obtaining the transfer function for a system repre- 
sented by eqixation (18) in response to a step input of the forcing func- 
tion is shown in appendix C. The transfer fimction is of the form 


Q 

8 


(p) 




+ J/ 




P + 


p'^ + bp + k 


(25) 


where a, b, k, and to are related by equations (20) and ( 21 ), and 
6 q is the magnitude of the step. 
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As a general rule, the exponential-approximation method. .of simu- 
lating transient data seems to offer the hest approach of-any of the 
methods reported in this paper when the response is a lightly damped 
oscillation to an approximate step. It is of- interest to note that 
this method may be used even thoi:igh the input is not approximated by 
exponential expressions, provided its analysis is restricted to the 
free-oscillation portions of the response. The coefficients k, b, 

J;l^ J2 eq.uations (18) and (19) may be obtained regardless of 

the form of input provided that they are obtained from a portion of-the 
time-response curve after the input has reached a steady-state value. 
This adaptation is pointed out by Shlnbrot in reference 16 and the 
method of application is reported therein. 

The foregoing method is useful only when the period and time to 
damp ta one-half amplitude may be read directly from the records . For 
highly damped systems where this is not possible and as an alternate to 
the foregoing approach, a least-sq.tiare 8 method for obtaining the period, 
damping, and other coefficients of equation (18) may be employed. 
Greenberg, in reference it-, disctisses extensively the application of the 
Prony method for fitting a sum of exponentials to a number of equally 
spaced ordinates. This method will, likewise obtain the transfer func- 
tion given by equation (25) . 

A meas\jre.of how closely the analytical expression represents the 
time response of the system during free oscillation _may be obtained by 
substituting the derived coefficients into equation (18) and allowing 
the time t‘ to vary. This substitution amounts to taking the inverse 
Laplace transform of the transfer function, once it has been determined, 
and returning the function Q(a)) to the time domain where it should be 
eqml to the original function of time q(t) . 

The time required for this method varies with the number of least- 
squares solutions required to obtain the four unknowns b, k, Jj_, 

and J2. The extreme case is the case where a least-squares solution 
is desired for all the coefficients of a response. In this instance if 
the input is considered to be a step and the response is described by 
2 k data points, the time required to obtain the transfer function of the 
system may be estimated at 8 to 10 man-hours of work and three separate 
least-squares solutions are required. Any reduction in the number of 
least-squares solutions will obvio^isly reduce the time required 
appreciably. 

Donegan-Peaxson method .- This method is appropriate for obtaining 
the transfer function from transient response to an arbitrary input, 
and, when only maniial computing techniques are available, offers a good 
degree of accuracy with a minimum of work. The method is presented in 



NA.CA TN 2997 


23 


reference 17 . In brief, a transfer function of one of the forms given 
by eq.uation (17) , for exan^ile. 


(Ap^ + Bp + C)q = 6(Ep + F) 


(26) 


is integrated twice and rearranged to give 



(26a) 

The equation is now in a form where the integrals may be calculated 
from q. and 6 which are known from a time history. The integral 
quantities in the foregoing equation are evaluated at some fixed time 
inteivals, starting with the initial control input, to form a series 

B C "P "P 

of simultaneous eq.mtions from which the — , — , and ^ coeffi- 

A A A A 

cients may be evaluated by the least-squares method. 


The expression of the transfer function in integral form is an 
Important point with regard to application of this technique in that 
the integration processes ane inherently more accurate than the differ- 
entiation process indicated in the normal form of the transfer functions. 
The integrals and the coefficients may be obtained by the matrix methods 
described in references I7 and I8. The use of higher-order terms in the 
numerator and denominator of the transfer function is possible but, in 
many cases, is iinnecessary because of the insignificance of their coef- 
ficients and is Impractical because of the large amoiint of additional 
work required. 

In the derivation of this method no restriction is made or Inplied 
that the forcing function (input) or transient response reach a steady 
state within the time limit considered. There appears to be, however, 
a practical limitation on the length of the record since, for any given 
short length of a curve, a large number of analytical expressions may 
be written that satisfy, with good precision, the conditions of the 
ciorve in the region considered. Obvio^^Bly, as the length of the record 
considered is extended, the expressions that adeq.uately define the curve 
become more limited until the correct expression is closely approached. 
In the analysis of eirbitrary inputs and responses that reached a steady 
state, this method produced excellent results over a large range of fre- 
q.uencies in a reasonably short length of time. When proficiency was 



2k 


NACA TN 2997 


obtained in the use of this method, including the matrix methods of 
integration and least squares, a congjlete frequency response required 
about 8 man-hours of work from raw data (averaging about 20 data points) 
to finished frequency- response curves (averaging 16 frequencies) . 

In the application of this method, a somewhat sin5)ler integrating 
matrix was iised that reduced the time of integration by about one -half 
over the method presented in reference 17 - The derivation of this 
integrating matrix by use of the relationships reported by Diederich in 
reference I8 is presented in appendix D. 

A check on the accuracy with which the time response is represented 
is afforded by the inverse Laplace transformation process described in 
reference 5 * A second method, which is suggested in reference 17, 
utilizes the eval\iated integrals of the output, the recorded input, and 
the transfer function. If the transfer function is to represent the 
system accurately, the response obtained by this check must-be equal to 
the original time response. 

For both the Donegan-Pearson and exponential-approximation methods, 
a second approximation to the determined transfer-function coefficients 
may be made by a procedure suggested by Shiribrot in reference 16. In 
cases where this refining procedure was attenpted, the process was 
lengthy and frequently did not afford better approximations because of 
failure of the method to converge. 


RESULTS AND C 0 I«>ARIS 0 NS OF METHODS 


The frequency responses as obtained from three aircraft of differ- 
ent types are used herein for' illustrative pvirposes._A summary of the 
mass and geometric parameters of these alrcraft^ogether with a sketch 
of their plan forms is presented in table I, as are the flight conditions 
for which the data were obtained. These aircraf't- will be referred to as 
a fighter, a transport, and a free-fall model. 

All methods of analysis were not applied to all three of these air- 
craft, but a comparison of the methods is made herein for the fighter at 
one flight condition. The conparison is made with the response in 
pitching velocity to an elevator step iioput; Time histories of the con- 
trol input and the response are sho'wn in figure 5 - The Donegan-Pearson, 
Prony, Fourier synthesizer, IBM, and manual Fourier methods were used to 
obtain the frequency response of the fighter from these time histories. 

A time Interval of 0.10 second was used for the manual methods, O.O5 sec- 
ond was used for the IBM method, and. O.O6 second was used for the Fourier 
synthesizer. The frequency-response curves thus obtained are shown in 
figure 6 , together with the faired cu 27 ve of figure 5 which was obtained 
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at the same flight conditions by the sinusoidal-response method. The 
frequency response of the fighter appears to he established to a gener- 
ally acceptable degree by any one of the several methods shown. The 
expression "acceptable degree" is expanded subsequently. 


Effect of Input Shape 

Inasmuch as the compatibility of results obtained by using the var- 
ious methods on a step input has been established, it is of interest to 
check the effect of this and other input shapes. 

A check of the effect of input shape on results obtained through 
use of the Donegan-Pearson method has been made. The step input illus- 
trated in figure 5 together with the approximately square and triangular 
inputs shown in figures 7 and 8, respectively, were lised for the inves- 
tigation and these inputs and their responses were analyzed at time 
intervals of 0.10 second. On each of these three figures, the accuracy 
with which the responses were represented by the transfer functions 
determined by the Donegan-Pearson method is shown by the data points on 
the response arrve. These points represent values of pitching velocity 
obtained by multiplying the integrated functions of equation (26a) by 
the derived transfer coefficients at the values of time indicated in 
figure 5. This procedure is described in more detail in reference I7. 

The frequency response of the fighter as determined by these trans- 
fer fiinctlons is shown in figure 9 compared with the faired curves of 
amplitude ratio and phase angle obtained from the sinusoidal -response 
method. The sln\isoldal-response method is included because it Involves 
a different test technique. These four sets of frequency-response cvirves 
appear to be in good agreement. Whether their agreement is to an 
"acceptable degree" may be illustrated by examining how closely they 
agree in the time domain when an identical control input is applied in 
each case. This process may be done manually through use of the inverse 
Laplace transformation; however, the Fourier synthesizer offers a machine 
method of obtaining the time response of a system described by the 
frequency-response curves to an approximate ramp or step . The Foiirier 
synthesizer was used in the present analysis and the inverse of the pro- 
cess described in the section entitled "Description and Disciission of 
Methods" was applied. 

The control input and time responses using the three frequency- 
response curves corresponding to the three input shapes investigated 
are shown in figure 10 . The curves show a maximum spread about I3 per- 
cent of steady-state value at steady state and a smaller percentage 
spread at peak overshoot. Thus, it appears that the determination of 
the transfer function is not particularly sensitive to the shape of the 
control input when the Donegan-Pearson method is used. 
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Effect of Harmonic Content of Input 

In the application of Fourier methods, the harmonic content of the 
input must he considered and has a predominant, effect on the results 
obtained. Harmonic content pertains to the relative magnitudes of the 
sine waves of various frequencies which make up the input or response 
shape and is essentially the amplitude of the Fourier transform of a 
function. 

In order to illustrate the hamnonic content of several inputs, fig- 
ure 11 shows the Fourier transforms of square, triangular, step, and 
impulse type of Inputs plotted against frequency. It can he seen that 
the harmonic content of— the square and triangular inputs go to zero at 
equally spaced increments of frequency, the spacing being dependent on 
the duration and shape of the input. For either- shape , doubling the 
duration of the input will halve the spacing between the frequencies of 
zero harmonic content. An error frequently encountered in the frequency 
domain when Fovirier methods eire used is caused by the harmonic content 
of the Fourier transform of the input closely approaching or reaching 
zero. When this condition occ\jrs, slight errors in the data cause the 
frequency-response curves to diverge and even to go to infinity at some 
frequency if the harmonic content of the input functions becomes zero 
at that frequency. 

An exanple of the distortion of the frequency-response c\arves due 
to low harmonic content was obtained in the analysis by Fourier methods 
of -the rectangular-pulse input and pitching- velocity response of the 
fighter as shown in figure 12. Because of- the length of the rectangular- 
pulse input used, the harmonic content of- both the input-and response 
closely approached zero at frequencies of about 8.5 radians per second. 
The discontinuity due to the lack of harmonic content is shown in the 
frequency- response plot~-of figure I 5 . An additional test utilizing an 
input which affords data having good harmonic content in the region of 
uncertainty wovild be desirable in order to Insure that-no legitimate 
secondeiry peak or other significant characteristic exist in that range 
of frequencies . 

When choosing inputs to be used in obtaining flight data for Fourier 
analysis, it is desirable to examine their harmonic amplitudes in light 
of expected instrument accuracies in order to select an input or series 
of inputs which will afford sufficiently accurate frequency-response 
data in the frequency range of interest. 

Since it is desirable to maintain large values of harmonic content 
over the entire frequency range, inputs approaching an impulse would 
appear most usable. In practice, however, control inputs of this type 
having suitable amplitudes must be maintained over a significant length 
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of time so that the airplane is disturbed sufficiently to insiare accurate 
measurement of the response. Thus, the transform of the resulting pulse 
vill often closely approach or reach zero at some frequency in the range 
over which the response is desired. 

The transform of the step input has the desirable feature of never 
becoming zero. !fe,vlng infinite amplitude at zero frequency, the trans- 
form decreases as the Inverse function of frequency and approaches zero 
amplitude as the frequency approaches infinity. In view of the rapid 
decrease in harmonic content with increase in frequency, however, it is 
sometimes difficult to maintain the accuracy of the frequency response 
to as high a value of frequency as desired. This effect may be seen in 
figure 11 by con^jaring the harmonic content of the triangular and step 
pulses from 1 to 8 radians per second. 

The trajisforms of several basic inputs together with the effect on 
harmonic content of distortion of these basic inputs are illustrated in 
reference 10. 


Effect of Record Length 

The advantage indicated for the so-called curve-fitting methods with 
regard to their ability to make a logical interpolation over frequency 
regions of low harmonic content woiild also appear to be applicable to 
extrapolation in either the frequency or the time domain. For example, 
the fact that the analytical form of the transfer function is assumed 
in advance for these curve-fitting methods would appear to afford possi- 
bilities for analyzing only a part of an input and response to establish 
the coefficients of the transfer fiinction whereas the Fourier analysis, 
by the nature of the limits of the Fourier integral, requires that a 
steady-state or a constant-amplitude oscillation be obtained. This 
apparent advantage of the curve-fitting methods in that the transient is 
not required to reach steady state has, in general, proven to have defi- 
nite practical limitations. In the application of both Fourier and curve- 
fitting methods, it has been found that time transients that do not closely 
attain steady state do not produce acciurate frequency-response data. 

Figure 12, which shows the time history of a rectangular-pulse ele- 
vator input and the response in pitching velocity of the fighter, may be 
used to illustrate these practical limitations . The data were analyzed 
in three stages. The response was first considered in the time interval 
from 0 ^ t ^ 0.70 second to be the response of the fighter to an approxi- 
mate step where the pitching velocity appeared to reach a steady state 
before the elevator was again disturbed. Time increments of O.O 5 second 
were used to obtain the frequency response of this portion of the time 
histories by the Prony, Donegan-Peanson, and manual Fourier methods. 
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(in this analysis, the Prony method required l6 man-hours; the Donegan- 
Pearson method required 9 man-hours.) In addition, check points were 
obtained by the manual Fourier method (by using Sin^ison's three-point 
rule of integration) with time intervals of 0.10, 0 . 025 ^ and 0.0125 sec- 
ond. The frequency-response results of this analysis are shown in fig- 
ure lif and, although all the methods closely agree, the frequency response 
appears to be quite different from that— indicated in figure 6 which was 
obtained for the same airplane at the same flight conditions. 

The source of this discrepancy was determined when the time histories 
shown in figure 12 were again analyzed by using the Donegan-Pearson method 
at 0 . 05 -second intervals to a time of l.ij-O seconds where the response 
still had hot approached a steady state too closely but the length of the 
record used had been doubled and the effective amplitude had been more 
than doubled. A third analysis was made by using the Donegan-Pearson 
method at 0.10-second intervals to 2.10 seconds at- which time its steady- 
state value was closely attained. The frequency response obtained by 
using each of the three record lengths is shown in figure I 5 together 
with the frequency response obtained for the fighter (by the Donegan- 
Pearson method) from figure 6. It may be seen that, when the first one - 
third of its response was analyzed, the record was short, and a steady 
state had not been reached; these factors precluded an adequately precise 
definition of the time response and an erroneous frequency response was 
obtained. When the length of the record was doubled, a more correct 
trend became apparent but, because a steady state had still not been 
defined, some fairly large discrepancies persisted, particularly with 
regard to the static value of the frequency response (the frequency- 
response curves of figure 6 being used as a basis for conqjarlson) . When 
the analysis included the entire response, even thovigh the time interval 
\ised in the analysis was doubled, a close agreement~with the frequency 
response obtained from the step input" was obtained. Reference I 7 rec- 
ommends that enough of the response time history should be taken to cover 
the natural period of the system. 


Other Causes and Effects of Errors 

In the determination of transfer functions from inputs and outputs 
having regions of low harmonic content, an advantage has been indicated 
to the approach of fitting an analytical expression to experimental data. 
In the authors' opinion, this curve-fitting technique, as compared to 
the Fourier- analysis, is of particoilar merit if there is reasonable con- 
fidence that the assumed analytical expression is of the correct form 
for the system being analyzed. In this manner, another condition (the 
form of the transfer function) is stipulated which the analysis must 
obey. In mathematical processes, the more conditions correctly stipu- 
lated, the more precise the results . On the other hand, errors in the 
transfer function or frequency response as obtained from the curve-fitting 
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methods due to either the wrong assumption of the form of the transfer 
fimction or due to the errors in the calculations are not readily appar- 
ent since the assumption of a given form will usually give variations 
that appear logical. However ^ as has been pointed out, certain checks, 
such as the use of the inverse Laplace transformation, are available for 
comparing the time response predicted from the transfer functions with 
the time-response curves from which the transfer functions were derived. 

In the Foinrier methods, inaccuracies are, in general, more readily 
discernible than in curve-fitting methods. In the use of Fourier methods, 
there has been found evidence of discrepancies attributable to three 
causes (as pointed out in ref. lO) : the lack of harmonic content of the 

Fo-urier integral, the use of too large time intervals in the time domain 
to afford accirracy in the frequency domain, and the Incorrect synchroni- 
zation of input and response data in the time domain. 

The first of these eirors has already been discussed in the consid- 
eration of the effect of input on Fourier methods and, as has been pointed 
out, is usually discernible by divergence of the curves in some small 
range of frequencies . 

The second of these errors, that of too large time intervals, is 
generally indicated by a scattering of the data points in the frequency 
domain where the magnitude of scatter ■osixally diverges rapidly with 
increasing frequency. Insight into the cause of this scatter may be 
seen in the characteristics of the Foinrier transform where, at'each fre- 
quency, the transient q.(t) is multiplied by a sine and cosine wave of 
unit aitplitude and where the resulting area under the two product curves 
determine the coefficients of the real and imaginary parts of the complex 
variable in the frequency domain. As frequencies greater than the natural 
frequency are Investigated, the differences in the positive and negative 
areas of the product curves grow smaller (compared with the magnitude of 
the individual areas) so that the effect of small errors is magnified. 
Thus, small Inaccuracies in the representation of the transient curve 
become more prominent as higher frequencies are Investigated and appear 
in the frequency domain as scatter. Several estimates of the frequency 
at which scatter will become inportant, for the different Fourier methods, 
based on the time Interval chosen, have been given in the section entitled 
"Description and Discussion of Ifethods." 

A typical occurrence of scatter due to the choice of too large a 
time lntei*val was obtained when the response of a free-fall model, the 
characteristics of which are given in table I, was analyzed at 0. 10-second 
intervals by the manual Fourier (numerical- Integration) method. The ele- 
vator input used and the response of the model, in angle of attack are 
shown in figure l6. The frequency response as determined by the numeri- 
cal manual Fourier, the Coradi harmonic analyzer, and the exponential- 
approximation methods of analysis are shown in figure 17 . The scatter 
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of points obtained by the manual Fotirier method of numerical integration 
occirrs at. frequencies greater than about 8 radians per second. Further 
analysis vith smaller time intervals of, say, O.O 5 second should provide 
better results in this region. 

In the study of missiles and free-fall models where low damping is 
generally encoxmtered, the use of the exponential-approximation method 
is pairticularly useful and requires a minimum of time. The response in 
figure 16 was analyzed by both the least-squares (Prony) method with 
0.10-second intervals and by direct coniputation. Both gave- identical 
coefficients and the frequency response obtained by using these coeffi- 
cients is also shown in figure 17 . The Donegan-Pearson method was 
attempted with this type of response but~did not produce coefficients 
that represented this lightly dairped system as exactly as it did for 
systems with high damping. The representation of the time response by 
the derived transfer function is illustrated in figure 16 where the 
inverse Laplace transform was applied to the transfer function obtained 
by the Donegan-Pearson and exponential-approximation methods to predict - 
the response to a step. The reason for this condition is that the trans- 
fer coefficients which primarily determine the period and damping of the 
oscillation are determined by the double integration and integration, 
respectively, of the output. The smoothing effect of these integration 
processes on any existing oscillation does not therefore enable accurate 
detection of the oscillation characteristics. 

Although scatter obtained by using the Fourier approach is indicative 
of inaccuracies, the converse does not apply Inasmuch as the absence of 
scatter in the use of Fourier- methods is not an indication of correctness. 
An illustration of this point was obtained in the analysis of the fre- 
quency response of the transport, tested under the conditfons given in 
table I. The response in pitching velocity to an elevator input is shown 
in figure I 8 . The manual (numerical- integration) Fourier method, analyzed 
at 0.20-second intervals, was used to determine initially the response 
at 1, 5^ 8 radians per second and these frequency-response 

points are indicated in figure 19 • Althoiagh the amplitude ratios and 
phase angles at the frequencies Investigated did not indicate scatter, 
when two additional frequencies ( 6.5 and 7.5 radians per second) were 
Investigated, the scatter became apparent” At a smaller time interval of 
0.10 second, the control input and time response were analyzed by using 
the Fovirier synthesizer, Donegan-Pearson, and again the manual (numerical 
integration) Fourier methods. Although small differences in the results 
are apparent in figure 19 ^ they do agree sufficiently well for most pur- 
poses. A check point at a frequency of 8 radians per second and with a 
0. 05-second time interval was made with the numerical-integration Fourier 
method. The resixLt essentially substantiated the value determined with 
0.10-second time intervals. 
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The third mentioned cause of error often incurred in the Fourier 
analysis was referred to as incorrect synchronization of the input and 
response data in the time domain. A shift in the correlation of the 
time scales between input and output, in turn, causes a change in the 
phase angles obtained in the frequency domain by an amo-unt directly pro- 
portional to the frequency. These erroneous values of lag or lead will 
be hard to detect regardless of whether Fourier or cvirve-fitting methods 
are used since incorrect but apparently logical frequency-response curves 
will usually occur. 

In order to avoid or reduce errors in the determination of transfer 
functions from flight data, it is highly desirable to use as large a 
control deflection as possible, but the magnitude of this control deflec- 
tion must also be compatible with the requirement that the stability 
parameters of the aiip>lane remain within their linear range. It also 
appears highly desirable to analyze responses from two or more input 
shapes at a given flight condition. A congjarison of the frequency 
response obtained from the same record by different methods has also 
proved a useftil check. 

As mentioned previously, examples indicating the times required to 
reduce flight data to frequency responses by the various methods dis- 
cussed herein are summarized in table II. The table enables the weighing 
of the time factor in choice of a method} however, the choice depends on 
other factors as well, such as availability of machine computing equip- 
ment and limitations inherent in the various methods as have been 
discussed. 


CONCLUDING REMUIECS 


In the foregoing study a number of considerations are indicated 
which pertain to the choice of methods in the determination of transfer 
fxmctions and frequency response from transient data. These considera- 
tions may be summarized as follows: 

In the methods which Involve the analysis of transient responses 
over short periods of time, a control input should be used that will 
afford (a) a close approach to a steady-state condition and (b) response 
amplitudes and harmonic content (covering the frequency range of interest) 
large eno\igh to give good instrument and reading accuracy yet small 
enough to keep the aircraft from departing from the flight condition for 
which the response data are desired. Vfhen flight data are analyzed, it 
appears highly desirable, as a check on the determined transfer function, 
to obtain responses from two or more input shapes at a given flight con- 
dition. A comparison of the frequency response obtained from the same 
record by different methods has also proved a useful check. 
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The method involving the oscillation of an aircraft through use of 
sinusoidal control Inputs requires a large amount of flight" time hut a 
relatively simple analysis. Satisfactory data may he obtained with a 
human pilot generating an approximate sinusoidal control input. For 
modes of aircraft motion which are nearly critically da,mped, the large 
amount of flight time can he reduced somewhat since continuous records 
may he taken while the frequency of oscillation is slowly changed to 
cover the frequency range desired. 

Two manual Fourier methods of analysing arbitrary inputs and their 
responses were investigated. In the first approach, analytical expres- 
sions within discrete intervals are fitted to the time response and 
input, and terms of the Fourier Integral are obtained analytically. In 
the second approach, the time response and input functions at selected 
times are first multiplied by the sine and cosine functions appearing 
in the Fourier Integral and the resulting product curves are integrated 
nimierlcally . The first approach appears to he basically more accurate 
when utilizing manml computing hut is not as flexible or as suited to 
machine calculations as the second approach. Special machines for 
acconplishing a Fo\arier analysis, such as the Fourier synthesizer and 
Coradi harmonic analyzer, afford a means for significantly reducing cal- 
culation time as compared to a manual approach. The two machines men- 
tioned give results conparahle to those obtained by the usual maniaal 
procedure in n-umerlcal Fourier analyses. Because of its principle of 
operation, the Coradi harmonic analyzer appears to be basically more 
accurate than the Fourier synthesizer and, in general, gave satisfactory 
results out to somewhat higher frequencies. The Coradi harmonic analyzer 
appears more capable of handling random variations than the Fourier 
synthesizer. The mechanical application of the nimierlcal Foiarier analysis 
through use of punch-card calculating machines (for exanple, IBM equip- 
ment) is a means for appreciably reducing calculation time. This approach 
appears to afford the possibility of obtaining greater precision in the 
calculations since the rapid computation makes feasible the use of smaller 
time intervals combined with more complicated and precise integrating 
formulas . 

The exponential-approximation and Donegan-Pearson methods establish 
an analytical expression for the transfer function which, in terms of the 
imaginary frequency variable, is continuous in frequency. The Fourier 
analysis, in contrast, does not furnish analytical expressions and gives 
values of frequency response only at selected frequencies . The Donegan- 
Pearson and Prony methods can be used satisfactorily when reasonable con- 
fidence exists as to the analytical forms of the transfer function (since 
the form must be assumed in advance). This approach will not, however, 
detect details of the frequency response that cannot“be approximated by 
the assimied form even if such characteristics exist in the time response. 
In contrast, Foiu?ier analysis will detect all details of the frequency 
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response ■which, are within the accuracy of the measurements and the calcu- 
lation procedure. The escponential-approximation method is best suited for 
lightly damped systems where the control input closely approximates a 
step or is of a form that possesses a simple Laplace transform. The 
longitudinal transfer f'unctions of oscillatory aircraft can often be 
determined by simple, direct computation from the meas'ured period, 
damping, steady-state value of the response, and phasing of the time 
response. Uonoscillatory transients that do not afford direct approxi- 
mation of the response may be approximated by a least-squares procedure 
known as the Prony method. The Donegan-Pearson method appears beat 
suited to systems that are not highly oscillatory and works well for 
inputs that are not necessarily represented by analytical expressions. 

When least-squares procedures must be applied in the Prony method, the 
Donegan-Pearson method generally affords shorter calculation time. When 
more than a few discrete frequencies are desired, the Donegan-Pearson 
method affords a more rapid approach than manual Fourier analysis. 

Fourier methods are more critical to the forms of the input than the 
Donegan-Pearson method and inputs should be chosen to a-roid regions of 
low harmonic content in the frequency range of interest. Although the 
Donegan-Pearson method appears to interpolate satisfactorily over regions 
of low harmonic content, it does not appear to be applicable to large 
extrapolation in either the frequency or time domain. In the use of 
this method, as in the Fourier methods, it is necessary to obtain data 
which closely approach the steady s-tate in order to predict accurately 
the low-frequency-response characteristics. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., J\me k, 1955* 
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APPENDIX A 

ESTIMTION OF THE ERRORS ENCOUNTERED USING A CONTROL INPUT 
THAT CHANGES IN FREQUENCY AT A CONSTANT RATE 


Much of the extensive flight-test time involved in obtaining 
freq.uency-response data hy using a sin\isoidal input can he eliminated 
if the sinusoidal input~is continuotisly changed in frequency at a slow 
rate. Since transients are constantly being introduced- and dying out 
because of this constantly changing frequency^ the error introduced by 
assuming that the response to this wave approximates the steady-state 
response to a constant-frequency wave may, at any given frequency, be 
a function of the natxiral frequency and daniping of the aiip)lane as well 
as the rate at which the frequency is changed. Investigation of the 
magnitudes of these errors in an^jlitude and phase "angle based on the 
response of the airplane to a constant-:^equency sinusoidal input was 
performed as follows: 

If a wave form of constant an^litude and constantly changing fre- 
quency is compared with a sinusoidal wave form having the same an^litude 
but constant frequency, there will occur, at a time herein assumed to be 
zero, a condition where the amplitude and instantaneous frequency of the 
two waves will be Identical. In the following derivation, the frequency 
at this instant is defined as cUq and the two waves are adjusted so that 

at this instant both waves are at their maximum amplitude. A second- 
order system is considered. The differential equation relating the 
response of the system to a constant frequency input~ls then 


2C \ 

— — + — D + l)x = cos co^t 


(Al) 


Similarly, for the varying cosine wave, this differential eqviation is 


^ d 2 ^ 25 




+ — ^ D + 1) X = cos(cDq + Ct)t 


(A 2 ) 


Gomparison of the two inputs show that their difference is effec- 
tively a phase difference which varies as the parabola 0 = Ct^. Since 
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a time-response solution of equation (A2) was too cumbersome to "be 
feasible, a linear -phase relationship was chosen that wo\ild approximate 
the paraholic-phase relationship and would afford a relatively simple 
time-response solution. A wave having a linear-phase difference with, 
a constant-frequency wave is, of course, another constant -frequency 
wave of a different frequency. Althoiigh a constant-frequency wave 
would not appear to be a good approximation to the variable -frequency 
wave under consideration, it will be shown that the difference between 
the original constant-frequency wave and the varying-frequency wave may 
be closely approximated by the difference between the two constant- 
frequency waves within the region of interest provided their frequencies 
are properly selected. 

The procedure \ised for establishing the frequency of the wave used 
in the approximation, in terms of the rate of change of frequency of 
the variable-frequency wave, is as follows: The actml phase-angle 

variation and the assumed approximation are illustrated in figure 20. 

The time interval over which the actual phase-angle variation was 
approximated was the interval which would enable transients Introduced 
by the varying-frequency wave to reduce to one-twentieth of their initial 

value ^”^1/20 = ^ parabola was approximated by a straight 

line chosen to pass through the parabola at the times t = 0 and 
2 

t = - — '^ 1 / 20 ’ approximation was chosen to balance the areas 

between the parabola and the straight line in the region of Interest. 
Substitution of the approximate phase angle into the varying-frequency 
input for the right-hand side of equation (A2) gives 


cos 



(A3) 


The relationships '^■]./20 ~ ~~ ^ Ct^) = 2C may be 


substituted into equation (A5) to give 
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For an example case where Dm = 1 radian per second per second and 
U)q = 8 radians per second, the original constant-frequency wave (eq. (A1) ) 

and the varying-frequency wave (eq. (A2)) are shown in figure 21 together 
with the wave used to approximate the phase-angle relationship between 
the original two (eq. (A4-) ) . In addition, the difference between the 
original constant-frequency wave and the varying-frequency wave is com- 
pared in fig\ire 21 with the difference between the original constant- 
frequency wave and the wave used in the approximation of the phase-angle 
relationship . — 

The errors in phase and anq)litude incurred in the response of the 
second-order system and caused by the use of a varylng-frequency-wave 
input instead of a constant-frequency-wave input may be obtained by 
determining the difference between the phase angles and amplitude ratios 
obtained from these two inputs . These errors are ^ven by 


cos 

cos 0)Qt 

si + Si D + 1 si + Si D + 1 


The substitution of D = JcOq for the response to cos cu^t and the 
substitution of D = jfov, - ^" 0 ^ "the response to cos /cu-. - SSL.jt 

V K/ - - V M 

will, at t = 0, give the error relationship in terms of amplitude and 
phase angle 
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Various values of dancing ratio freq.uency ratio ^^<3. the 

rate of change of frequency ratio were substituted into the 

foregoing relations in order to obtain plots of phase error and ampli- 
tude error (related to the response to a pure sinusoidal input) at zero 
time. Figure 22 presents plots of these errors over a range of values 
of parameters pertinent to most aircraft- This figure indicates that 
the errors Increase rapidly with a decrease in the value of damping 
ratio below O. 707 . At low values of damping ratio, excessively large 
errors will be obtained unless the frequency is varied at an extremely 

low rate in the vicinity of = 0.01. The greatest errors in all cases 

appear to occur in the vicinity of the natural frequency, the errors 
approaching zero at high and low values of the frequency ratio. 

For the tests presented in the body of this paper, the airplane 
tested had a damping ratio of about 0-7- The rate of change of fre- 
quency for these tests was not constant but rather was held constant at 
one frequency for several oscillations before progressing to a new fre- 
quency. However, averaging over a range of frequencies gave an average 

value of of O.O 6 . Figure 22 indicates an error in an^plitude 

and in phase angle of less than 10 percent for the airplane tested at 
these conditions. This error falls within the scatter shown in figure 5* 
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APPEKDIX B 

METHODS OF FOURIER ANALYSIS WHEREIN THE TRANSIENTS 
ARE REPRESENTED BY A SERIES OF POLYNOMIALS 


A representation of an approximate step-control input and tiie time 
response of pitching velocity of the fighter as shovm in figure 5 involved 
the division of the input^nto three intervals of time and the division 
of the response into five intervals of time. 

The input, which reaches a steady-state value of- 0.74° at 
t = 0.10 second, was represented within the intervals hy the equations 

5Q(t) = 20 t^ (0 ^ t ^ 0 . 05 ) 

= 13 . 8 t - 0.64 ( 0.05 $ t ^ 0 . 10 ) 

62 (t) = 0.74 ( 0.10 < t S «) 


The Fourier transform of the input- is then 


5(03) 



.2 -Jest,, 
t e dt 




0.10 

0.05 


(I5.8t - 0.64)e"'^^dt + 


0.74 r e-^\t 
Jq.iq 


13 . - ll.Se-Q-Q^J^ ^ j 40(1 - 
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or in trigonometric form 


5(cd) = cos O.lCtD - cos O.O 5 CO - ~ sin 0.05tu + 

Oi^ 


sin O.lOoii + — (1 - cos 0.05cu) (B1) 


0.2 

0.2 

sin 0.05oi 

13.8 


or 


or 


where the substitution of selected values of freq.uency oj -will afford 
the real and Imaginary coefficients of the Fourier transform at each 
frequency chosen. 

The time response of pitching velocity, shown in figure was 
divided into discrete intervals and the coefficients of the cubic 
equation 


D9(t) = at^ + ht^ + ct + d 


were fo\and as follows; 

For the time interval 0 ^ t S 0.10, inspection of the curve indi- 
cates that it may he closely approximated hy a cubic without lower-order 
terms so that b .= c = d = 0 and, at t = 0.10 second, the relationship 
is written 


D9(t) = at^ 


0 . 0014-9 = a(O.lO) 


5 


k.9 = a 


The equation for this time interval becomes 
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For the time interval 0.10 t < 0 . 50 ^ a q^uadratic representation, 

since the transient in this interval does not indicate the need of a 
cubic representation, will be assumed where the coefficients are found 
by the solution of the following equations: 


D 9 (t 

De(t 

D 9 (t 


0.10) 

0.0049 = (O.lO)^b + 

(O.lO)c + 

0 . 20 ) 

0.0489 = ( 0 . 20 )^ + 

( 0 . 20 )c + 

0.50) 

0.0735 = (0.50)^b + 

(0.50)c + 


d 

d 

d 


The solution by "successive subtraction" is illustrated here since equal 
time Intervals were used. 


O.OOli -9 - O.OI189 = (0.01 - 0 . 04 )b + (0.10 - 0.20)c + 
0.0489 - 0.0755 = ( 0.04 - 0.09)b + (0.20 - 0.50)c 
0.0440 = 0.05b + 0 . 10 c 
0.0244 = 0,05b + 0.10c 

0.0440 - 0.0244 = (0.05 - 0.05)b + (O^l^-'-Ovlliic^^ 


0.0196 = -0.02b 
-0.98 = b 

0.754 = c 


-0.0587 = d 
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The equation for this time interval becomes 

D03_(t) = - 0.0587 

For the time interval O.3O < t < O.8O, the slope of the transient 

is zero at a time of 0.55 second. Use of this condition is desirable 
in evaluating the constants over this interval; therefore, 

D 9 (t = 0.30) 0.0733 = (0.30)^a + (0.30)^ + (0.30)c + d 

De(t = 0.55) 0.0856 = (o.55)5a + (o.55)^b + (0.55)c + d 

D 0 (t = 0.80) 0.0793 = ( 0 . 80 )\ + ( 0 . 80 )^b + ( 0 . 80 )c + d 

D^ 0 (t = 0.55) 0 = 3 ( 0 . 55 ) 2 a + 2 ( 0 . 55 )b + c 

and the equation for this time interval becomes 

U02(t) = 0 .l 92 t^ - O.J4^55t^ + 0 . 3379 t + 0. 00862 

For the time interval O.8O ^ t ^ 1 . 40 , the solution of 

D 0 (t = 0.80) 0.0795 = ( 0 . 80 )^a + ( 0 . 80 )% + ( 0 . 80 )c + d 

D 0 (t = 1 . 00 ) 0.0725 = (l. 00 ) 5 a + (1.00)% + (1.00) c + d 

D 0 (t = 1 . 20 ) 0.0665 = (l. 20 )\ + (1.20)% + (l. 20 )c + d 

D0(t = 1.14-0) 0.0636 = (l.l4-0)5a + (l. 40 )% -4- (l. 40 )o + d 

gives an equation for this time interval oi 

D05(t) = 0.035l4-t^ - 0 , 09125 t^ + O.Ol4-25t + O.O8565 
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For the time Interval 1.40 ^ t ^ k>, the equation for a constant 
value from steady state to infinity becomes 

D9(t = 1.40) 0.0636 = d 


Therefore , 


D0ij^(t) = 0.0656 


The Fourier transform of the response is now evaluated by using the 
foregoing analytical expressions by s-ummlng the following integrals: 


D9(cu) = r D0(t)e *^^dt 
'^0 


pO. 10 , . ro.30 

= 4.9 / t^e"*^ dt - / (0.98t^ - 0.754t-+ 0.0587)e"'^’*^dt + 

^0 "^ 0.10 



(0.l92t5 - 0.4655t2 + 0.5579t + 0 . 00862 ) e'^^^dt +- 



(0.0354t^ - 0 . 09125 t^ + 0.0425tr-+ 0. 08565 ) e"'^\t + 


0.0636 r e"'^’*^dt 
^1.40 
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The integration of these tenas (given in text after eq.. (6)) Leads to 
the following relations: 


D9(cu) 


^ 8.8 _ p 8.8 

03^ 03^ V 03^ 



(cos O.lCtD - J sin O.lCtD) + 






(cos O.JCto - j sin 0.5Cto) - 


i.88 


0 . 003 ]^ - (cos 0.8013 - j sinO. 8013 ) - 


0.4248 


.004-^ + (cos 1.40 i 3 - j sin 1.40i3) (B2) 


where the substitution of selected values of 03 affords the real and 
imaginary coefficients of the Fourier transfomi of the pitching- velocity 
response of the fighter. The relationship between the output and the 
input is given by eqiiations (4), (6), and ( 7 ) and a plot of these rela- 
tionships at several values of frequency are shown in figure 6. 
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APPENDIX C 

APPLICATION OF THE JLAPLACE TE?ANSFQRM TO A RESPONSE EQUATION 


The response equation for the normal acceleration, 
Prony method described by Greenberg in reference 4, is 


as used in. the 


n(t) = ngg + sin cut + Jg cos £tot) (Cl) 

where rigg is the normal acceleration at steady state. The equation for 
a step control input of magnitude 6 q is as follows: 


5(t) =60 (t ^ 0) (02) 


The Laplace transform, indicated by the operator p when applied 
to equation (Cl) , gives 


n(p) 


ngg J 2 (p - a) 

P (p - a)^ + u>^ (p - a)^ + CD^ 


(C5) 


and equation (C2) becomes 

sCp) = . (oit) 

P 

By definition, the transfer function is the ratio of the Laplace transform 
of the output to the Laplace transform of the input _( Initial conditions 
zero) 


n(p) 

6(p) ^o 


+ JgPCp - a) 
4 * — 

p^ - 2ap + a^ + 


8 ^ 


i^ss + 


J2P^ + - J2a)p 


p^ -h bp + k 


(C5) 
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where the suhstitutlon of b = -2a nnfl k = has been 

equation may be rearranged to agree with the form of equations (17) 
that 


2 ^ss + ^2 


|(P) = 


+ P 


bngg + - J2a kn. 


ss 


+ bp + k 


iv5 


The 

so 


(C6) 


where the substitution of jco for the operator p will produce the 
frequency-response relationships. 
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APPENDIX D 

AN INTEGRATING MATRIX 


The Donegan-Pearson method (ref. I7) of analysis of transient 
responses suggests a matrix solution which, if followed, requires some 
knowledge of integrating matrices. An integrating matrix believed to 
be somewhat easier to use than the one su gg ested in reference IJ is 
presented here together with its derivation. 

If an arbitrary time-curve is chosen and divided into equal inter- 
vals of time, then by Simpson's rule a parabola may be described through 
three adjoining points _ 



By the use of equation (a 4 ) of reference I8, 
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Solution of equation (D1) for values of the Integers n = 0, 1, 2, 5, 
4, . . . m gives 



y dt = ^(5yn_4 + 157x1-5 ^^n-2 + 12yn-l + ynn) 

i^n-il 


f ^ y dt = ^(5yn_ni + 15yn-(ni-l) + '^n-im-2) + • • • + 

■^n-l 

12yn-l + 77x1 - yxH-l) 

Written in tabular form (-without fixing a value for At ) , the coef- 
ficients of y appear in the integrating matrix in the form 


t 

0 

At 

2 At 

— 
5 At 

4 At 

5 At 

6 At 

7 At 

0 

0 

0 

0 

0 

0 

0 

0 

0 

At 

O.IU6667 At 

0.666667 At 

-0.0^333 At 

0 

0 

0 

0 

0 

2 At 

A16667 At 

1.083335 At 

.583555 At 

-0.083333 At 

0 

0 

0 

0 

5 At 

.IU6667 At 

1 . 083335 - At 


.583355 


0 

0 

0 

4 - At 

A16667 At 

1.083333 At 

1.000000 At 

1.000000 At 

.585555 At 


0 

0 

5 At 

.lil 6667 At 

1.083335 At 

1.000000 At 

1,000000 At 

1.000000 At 

.585555 At 

-0.083333 At 


6 At 

.416667 At 

1.083555 At 

1.000000 At 

1.000000 At 

1.000000 At 

1.000000 At 

.583555 At 
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The use of this integrating matrix sln^ily requires the accimrulative 
summation of— “ 


(0.416666 + (0.666666 At)yj^ - (0.085555 At)y^^ 


which gives the area lying between n - 1 
already foiand from 0 to n - 1. 


n added to the area 


and 
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TABLE I 


CHAEACTERISTICS AITO FLIGHT CONDITIONS OF AIRCRAFT 
USED FOR LONGITUDINAL TESTS 








TABLE n 


VJl 

M 


EXAMPLES OF TIMES KB(3aiRED TO OBTAIN FREQUENCT DATA 
FROM TRAHSIMT DATA BY VARIOUS METHODS 


Method of analysis 

Number of 
functions 
analyzed 

Number of 
data points 
per function 



Number of 
frequencies 
analyzed 
per function 

— 
Time required 
to obtain 
transfer 
functions, hr 

Time required 
to obtain 
frequency 
response, hr 

Sinusoidal 

2 

— 

16 


2 

Fourier : 






Analytical integration 

2 

15 

11 


12 

Numerical Integration 

2 

2k 

8 


12 

Punch- card coniputer (IHM) 

2 

241 

12 


10 

Fourier synthesizer 

2 

24 

®12 


4 to 8 

Coradi harmonic analyzer 

2 

Continuoua 

10 


8 



curve 




Curve fitting: 






Emponential approximation 






Direct confutation 

^1 

. 24 

® 15 . 

0.5 

2-5 

Prony 

^1 

24 

^15 

6 to 8 

8 to 10 

Donegan-Pears on 

2 

20 

H6 

6 

8 


^'Results were obtained In form continuoua with frequency hut numerical evaluation was 
made at number of freq.uencies shown. 

^^thod Is primarily applicable to in^julse or step ir?>uts; therefore^ a n alysis of the 
input function is not required. 
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Time, t, sec 


Figure 2.- Film records of sinusoidal elevator input and response of 
filter in pitching velocity at M =s 0.6 and = 10,000 feet. 
Average an5)litude, ^ (Ap + A2 + A5); average period, i (Pi + P2 + P5); 

average lag, i (ti + Tg + 13) . 
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Figure 3-- Sinusoidal-response data points of fighter airplane at 
M = 0.6 and hp = 10_,000 feet. 




Nornuil acceleration, n Elevator Input, 0 ^, deg 



Figure 4.- Relatlonahip of leg t to phase angle 0 in a time response 


to a step elevator Input. 


y) = at ^ = ac(i 
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Pitching velocity, DO, radians/sec Elevator in 


to 
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'O 
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a 


gure 5*- • 
response 
points i] 
obtained 









Phase angle, j6, deg Amplitude ratio 



Fre 

6.- Frequency respon 
ilevator deflection a 
everal methods. 




Pitching velocity, D0, radians/eec Elevator input, 6e, deg 
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Time, t, sec 


Figinre 7.- Time history of elevator square-pvO.se input and pitching 
velocity response of fighter at M = 0.6 and hp = 10,000 feet. 
Circled points indicate the response calculated frcan transfer 
coefficients obtained from Donegan-Pearson method. 



Pitching velocity, D0, radians/sec Elevator input, 6®, deg 


6o 


NACA TN 2997 




Figure 8.- Time history of elevator triangular -pulse input and pitching- 
velocity response of filter at M = 0.6 and hp = 10,000 feet. 

Circled points Indicate the response calculated from transfer coeffi- 
cients obtained from Donegan- Pearson method. 




0 U Q 

Frequency, w, r 


Conrparisoii of frequency-response 
and hp = 10,000 feet as ottaini 

g tlie Donegan-Pearson and sinusoii 





Time, t, sec 


Figure 10.- S imu lated elevator Input and time respons 
pitching velocity as defined hy several frequency- 





An^Jlitude or Fourier trarsfonn as a rimction of A 
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Time, t, seo 


Figure 18.- Time history of elevator rectangular-pulBe input and pitching- 
velocity response of fighter at M = 0.6 and h^, = 10,000 feet. 
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Figure 16.- Time history of elevator step input and. angle- of-attack 

response of the free-fall model at M = O.725 and hp = 52,000 feet. 
The responses, as computed hy the Donegan-Pearson and exponential- 
approximation methods, are Indicated. 
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phase angle 



Figure 20 .- Straigjht-line approximation to the paraholic-phaBe- angle 
variation vtiich exists "between a constant-frequency wave and a wave 
whose frequency changes at a constant rate. 
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Figure 21.- Time history of cosiae wave forms in the vicinity of zero 
time for a case where Du) = 1 reidian per second, per second and 
cdq =; 8 radians per second. 




